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In recent years, advances in computational power and spatial data analysis (GIS, 
remote sensing, etc) have led to an increase in attempts to model the spread and 
behaviour of wildland fires across the landscape. This series of review papers endeav- 
ours to critically and comprehensively review all types of surface fire spread models 
developed since 1990. This paper reviews models of an empirical or quasi-empirical 
nature. These models are based solely on the statistical analysis of experimentally 
| obtained data with or without some physical framework for the basis of the re- 

lations. Other papers in the series review models of a physical or quasi-physical 
nature, and mathematical analogues and simulation models. The main relations 
of empirical models are that of wind speed and fuel moisture content with rate 
of forward spread. Comparisons are made of the different functional relationships 
<^ . selected by various authors for these variables. 
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Abstract 



History 

An empirical model is one that is based upon observation and experiment and not on 
theory. Empiricism has formed the basis for much of the scientific and technological 
advances in recent centuries and generally provides the benchmark against which theory 
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is tested. The study of fire and combustion in general was mainly an empirical endeavour, 
directed primarily toward application of combustion to industrial processes (for example, 
the industrial revolution of the late 1700s-1800s), until early in the previous century when 
the physical or theoretical approach had matured to the point of providing significant 
advances in understanding and prediction. The development of physical understanding of 
other forms of combustion (i.e. unintentional or uncontrolled fire) in general and wildland 
fires in particul ar, did not occur, however, until only very recently (in the last few decades) 
flSullivanl . [2007h . 



While there had always been a great general interest in unintentional fire in urban settings 



([Williams! Il982l ). for instance, the Great Fire of London in 1666, or the Chicago Fire 



on October 1871-prevention, control, prediction-unintentional fire in wildlands received 
much less attention, mainly due to the relatively little impact such fires have on the general 
populace. The study of the behaviour of fires in wildland regions has traditionally been 
driven by the needs of those practitioners involved in wildland resource management- 
foresters for the most part-for whom understanding this natural phenomenon was critical 
to the success of their work. 

Despite the fact that practically no region of the world (except for Antarctica) is free from 
such fires, much of the work in this fiel d was galvan ised in the United Stat es follow i ng the 
dev astating 1910 fire s in th e mid- west (jPynd . 120011 ). where workers such as iHawleyl (119261 ) 
and iGisbornd (119271 . Il929l ) pioneered the notion that understanding of the phenomenon 
of wildland fire and the prediction of the danger posed by a fire could be gained through 
measurement and o bservation and theoreti c al co nsider ations of the factors that might 
influence such fires. Curry and Fons ( 1938 . 1940l ). and Fons ( 1946 ) brought a rigorous 
physical approach to the measurement and modelling of the behaviour of wildland fires 
that set the benchmark for wildland fire research for decades following. 

In addition to the work conducted in the US, through the Federal US Forest Service and 
State agencies, other countries became increasingly involved in wildland fire research, 
primarily through their forest services-the Canadian Forest Service, the Commonwealth 
Forestry and Timber Bureau (later absorbed into the Commonwealth Scientific and In- 
dustrial Research Organisation (CSIRO) in conjunction with various state authorities in 
Australia-although many other countries such as South Africa, Spain, Russia, France, 
Portugal to name a few, have also had significant impact on wildland fire research. 

Since the early 1990s, European Union countries have committed significant funds towards 
wildland fire research, resulting in a boom period for this research in mainly Mediterranean 
countries and a major shift in focus away from the pioneering three (US, Canada and 
Australia). 

During the past two decades, the direction of much of the wildland fire research has been 
toward the use of fire as a resource managemen t tool in the form of hazard reduction burn- 
i ng or t he study of ecological effects of fire (e.g. iGill et al.l (119811 ) ; iGoldammer and Jenkins 
(Il990h : lAbbot and Burro"w^ (bOQSl ). 



Empirical modelling 



The focus of empirical modelling of wildland fire in the past has been on the determination 
of the key characteristics used to describe the behaviour of the fire. These generally have 
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been the rate of forward spread (ROS) of the head fire (that portion of the fire perimeter 
being blown downwind and normally of much greater intensity that the rest of the fire 
perimeter), the height of the flames, the angle of the flames, and the depth of flames at 
the head, although other characteristics such as rate of perimeter or area increase may 
also be of some interest. 

While observations of wildfires or fires lit intentionally for other purposes (such as hazard 
reduction or prescribed fires) have been used in the development of empirical models of 
fire behaviour, the predominant method has been the lighting of 'experimental' fires- 
fires whose only purpose is that of an experimental nature. This method can be divided 
into four parts. Firstly, the characterisation and quantification of the fuel and terrain 
in which the fire will lit (the slowly varying variables, which has included fuel load, fuel 
height, moisture content, bulk density, combustion characteristics, slope, etc.). Secondly, 
the observation and measurement of the atmospheric environment (the quickly varying 
variables, wind speed and direction, air temperature, relative humidity, etc.). Thirdly, the 
lighting, observation and measurement of the fire itself (its speed, spread, flame geometry, 
combustion rate, combustion residues, smoke, etc.). Fourthly, the statistical correlation 
between any and all of the measured quantities in order to produce the model of fire 
behaviour. Many workers have chosen to limit or control the possible natural variation in 
many quantities by conducting experimental fires in laboratory conditions which aids in 
the analysis of such fires. 

The primary use of such models has been to estimate the likely spread in the direction of 
the wind (and potential for danger to firefighter safety) for suppression planning purposes, 
much of which has traditionally been conducted in the form of simple 'back of the envelope' 
calculations for plotting on a wall map. Due to this simple need, empirical fire spread 
models have traditionally been one dimensional models in which the independent variable 
that is predicted is the rate of forward spread of the head of the fire in the direction of 
the wind. The rather pragmatic nature of these models, their relatively straightforward 
implementation, their direct relation to the behaviour of real fires, and, perhaps most 
importantly, their development by for the most part by forestry agencies for their own 
immediate use, have meant that empirical fire spread models have gained acceptance with 
wildland fire authorities around the world and to varying degrees form the basis for all 
operational fire behaviour models in use today. 



Operational models 



In the United States, the quasi-empiric al model of iRot 



National Fire Danger Rating Syste m (|Deeming et al. 



lermgj (Jl972h forms the basis of the 



19771 ; iBurganl . 119881 ) and the fire 



behaviour prediction tool B EHAVE (Andrews . 19861 ). This model is based on a heat bal- 



ance model first proposed by lFransdenl (I197ll ) and utilised data obtained from wind tunnel 
experiments in artificial fuel beds of varying characteristics and from Australian field ex- 
periments of grassfires in a range of wind speed condi t ions t o correlate fire behaviour with 
measured input variables. The model of Rothermel ( 1972 ) and associated systems have 
been introduced to a number of countries, particularly Mediterranean Europe. 

In Australia, the pred ominant operational fi re spread pre diction sys t ems h ave been the 
Mc Arthur Grassland (jMcArthurl . 1 19651 . 119661 ) and Forest (jMcArthurl . 119671 ) Fire Danger 
Rating Systems (FDRS), and the Forest Fire Behaviour Tables for Western Australia 
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(comm only called the Red Book) (jSneeuwjagt and Peetl . Il985h . based on the work of IPeet 
(119651 ). Both McArthur's systems and the Red Book are purely empirical correlations of 
observed fire behaviour and measured fuel and environmental variables from mainly field 
experimental fires augmented by well-d ocumented wi l dfires . More recently , the CSIRO 
Grassland Fire Spread Meter (GSFMj fjCSIROl. Il997l : IChenev and Sullivanl . Il997h based 
on the empirical modelling of Cheney et al. f 19981 ) has replaced the McArthur Grassland 
FDRS as the preferred tool for predicting fire behaviour in grasslands. This, too, is based 
on field experimentation and documented wildfire observations. 



In Ca nada, the quasi-empirical Fire Behaviour Prediction (FBP) System ( [Forestry Canada Fire Danger Grc 
19921 ) forms part of the Canadian Forest Fire Danger Rating System (CFFD RS) and is the 
culm i nation of 60 years of research effort in fuel moisture and fire behaviour (IVan Wagner , 
19981 ; iTay lor and Alexander! . 120061 ). Almost 500 fires were used in the construction of the 
FBP system, of which approximately 400 were field experiments, the remainder well- 
documented observations of prescribed and wild fires. The CFFDRS has been introduced 
and implemented in a number of countries, including New Zealand, Mexico and several 
countries of south-east Asia. 



The main characteristic of all but the CSIRO GSFM is that these systems were based 
primarily on small (<1 ha) experimental or laboratory fires and augmented with wild- 
fire observation^. T he series of experiments upon which the CSIRO GFSM was based 
( Cheney et all Il993l ) was the first to use experimental burning plots of which the smallest 
was 1 ha (See review of this model below). 



Background 



This series of review papers endeavours to comprehensively and critically review the ex- 
tensive range of modelling work that has been conducted in recent years. The range 
of methods that have been undertake n over the years represents a continuous spectrum 
of possible modelling (jKarplud . 119771 ). ranging from the purely physical (those that are 
based on fundamental understanding of the physics and chemistry involved in the be- 
haviour of a wildland fire) through to the purely empirical (those that have been based on 
phenomenological description or statistical regression of fire behaviour). In between is a 



continuous meld of approaches from one end of the spectrum or the other. IWeberl (Il99ll ). 
in his comprehensive review of physical wildland fire modelling, proposed a system by 
which models were described as physical, empirical or statistical, depending on whether 
they accounted for different modes of heat transfer, made no disti n ction between different 
heat transfer modes, or involved no physics at all. IPastor et al.l (120031 1 proposed model 
descriptions of theoretical, empirical and semi-empirical, again depending on whether the 
model was based on purely physica l unders t andin g, of a statistical nature with no physical 
understanding, or a combination. iGrishinl (119971 ) divided models into two classes, deter- 
ministic or stochastic-statistical. However, these schemes are rather limited given the 
combination of possible approaches, and, given that describing a model as semi-empirical 
or semi-physical is a 'glass half-full or half-empty' subjective issue, a more comprehensive 
and complete convention was required. 



3 Another interesting characteristic is that the CSIRO GFSM model was the only one published in a 
peer-reviewed journal; all the others were published as technical reports by the associated organisations. 
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Thus, this review series is divided into three broad categories: Physical and quasi-physical 
models; Empirical and quasi-empirical models; and Simulation and Mathematical analo- 
gous models. In this context, a physical model is one that attempts to represent both the 
physics and chemistry of fire spread; a quasi-physical model attempts to represent only 
the physics. An empirical model is one that contains no physical basis at all (generally 
only statistical in nature), a quasi-empirical model is one that uses some form of physical 
framework upon which to base the statistical modelling chosen. Empirical models are fur- 
ther subdivided into field-based and laboratory-based. Simulation models are those that 
implement the preceding types of models in a simulation rather than modelling context. 
Mathematical analogous models are those that utilise a mathematical precept rather than 
a physical one for the modelling of the spread of wildland fire. 

Since 1990 there has been rapid development in the field of spatial data analysis, e.g. 
geographic information systems and remote sensing. Following this, and the fa ct that 



there has not been a comprehensive critical review of fire behaviour modelling since [Weber 



(1199 ll ). I have limited this review to works published since 1990. However, as much of 
the work that will be discussed derives or continues from work carried out prior to 1990, 
such work will be included much less comprehensively in order to provide context. 



Previous reviews 



Many of the reviews that have been published in recent years have been for audiences other 
than wildland fire researchers and conducted by people without an established background 
in the field. Indeed, many of the reviews read like purchase notes by people shopping 



Perrv 


( 


1998); 


Pastor et al. ( 



endeavouring to be comp rehensive, hav e offer ed only superficial and cursory inspections 
of the models presented. Morvan et al.l ( 2004 ) take a different line by analysing a much 



broader spectrum of models in some detail and conclude that no single approach is going 
to be suitable for all purposes. 

While the recent reviews provide an overview of the models and approaches that have been 
undertaken around the world, mention must be made of significant reviews published much 
earlier that discussed t he pr ocesses in wildland fire propagation themselves. Foremost is 
the work of IWilliamsl (119821 ) which comprehensively covers the phenomenology of both 
wildland and urban fire, the physics and chemis try of combustion, and is re commended 
reading for the beginner. The earlier work of Emmons ( 1963 . 1966 ) and Lee ( l~9~7~2h 
provid es a sound background on the advances made during the post-war era. iGrishin 
(119971 ) pr ovides an extens i ve rev iew of the work c ondu cted in Russia in the 1970s, 80s 



and 90s. I Chandler et al.l (119831 ) and iPyne et al.l (119961 ) provide a useful review of the 



forestry approach to wildland fire research, understanding and practice. 

The first paper in this series discussed those models based upon the fundamental principles 
of the physics and chemistry of wildland fire behaviour. This particular paper will discuss 
those models based directly upon only statistical analysis of fire behaviour observations 
or models that utilise some form of physical framework upon which the statistical analysis 
of observations have been based. In this paper, particular distinction is made between 
observations of the behaviour of fires in the strictly controlled and artificial conditions of 
the laboratory and those observed in the field under more naturally occurring conditions. 
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The last paper in the series will focus upon models concerned only with the simulation of 
fire spread over the landscape and models that utilise mathematical conceits analogous 
to fire spread but which have no real-world connection to fire. 



Empirical models 

The following sections identify and discuss those empirical and quasi-empirical surface- 
only fire spread model s that app e ared in the literature since 1990. It is interesting to 
note the observation of ICatchpok (l2000h that the majority of new models that have been 



developed in recent years have been the result of efforts to initially develop and validate 
local fuel models required for the implementation of the BEHAVE (based on Rothermel) 
fire behaviour prediction system. Many researchers obviously felt that it was far easier to 
start from scratch with a purpose built model than to try to retrofit their local conditions 
into an existing model. Table [1] summarises the empirical models discussed in this review. 

Due to the varied nature of the empirical models presented here, including the fuels 
and weather conditions under which the data for the construction of the models were 
collected, the size and number of experimental fires and purposes for which the models 
were developed, it is difficult to compare them side by side. One possible method is 
the relationship between rate of forward spread (ROS) and wind speed. Wind speed is 
widely accepted as being the dominant variable determining the forward speed of a fire 
front. The reasons for this are cause for significant debate, ranging from the reduction in 
angle of separation of flame to unburnt fuel to increased turbulent mixing of combustants. 
Regardless of the mechanics of the process, the empirical approach to modelling fire spread 
must cater for this process and is manifested in the functional form chosen to represent 
it. Fuel moisture content (FMC) is also a key variable in determining rate of spread and 
this is also discussed. Fire spread models for fuel layers other than surface fuels, such as 
crown fires or ground fires, are not covered. 



Canadian Forest Service (CFS) - Acceleration (1991) 



While the Canadian Forest Service (CFS-accel) work (IMcAlpine and Wakimotol . 119911 ) 
is not a model of fire spread as such, it does address a major concern of fire spread, 
namely the acceleration in rate of fire spread from initiation. The assumption is that a 
fire will attain an equilibrium rate of spread for the prevailing conditions (the prediction 
of which is the primary aim of all fire spread prediction systems discussed here). The 
form of the function for the time to reach th i s equ ilibri um ROS is assumed to be expo- 



nential based on models propose d by Cheney ( 1981 ) and Van Wagner f 19851 ) (as cited by 



McAlpine and Wakimotc 



p ropose d 

3 JiaaS)). 



29 experimental fires were conducted in a wind tunnel with a fuel bed 6.15 m long by 0.915 
m wide consisting of Pinus ponderosa needles or excelsioi0 of varying fuel load and bulk 
density. Four wind speeds (0, 0.44, 1.33 and 2.22 m s _1 ) measured at mid-flame height 
were used. Temperature and relative humidity were held constant at 26.7°C and 80%. 



4 Excelsior is wood shavings cut into long thin strands 
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Equilibrium ROS was assumed to occur after 2.0 m of forward spread and determined 
using linear regression of averaged fire location and time measurements. 

Acceleration was modelled as an allometric (power law) function asymptoting to the equi- 
librium ROS with two coefficients, one based on the equilibrium ROS (which eliminated 
differences in fuel properties and integrated all other burning condition variables) and the 
other on the wind speed. The model was found to well represent the laboratory data 
but observations of elapsed time to equilibrium ROS did not coincide with point source 
field observations of other authors which were much greater and also dependent on wind 
speed. 



CALM Spinifex (1991) 



The Western Australia Department of Conservation and Land Management (CALM) 



Spinifex model (IBurrows et al.l . Il99ll ) was developed from 41 experimental fires conducted 



in predominantly spinifex (Triodia basedowii and Plectrachne schinzii) fuels on relatively 
flat sand plains. These fires were lit using drip torches to create lines of fire up to 200 m 
long perpendicular to the wind direction. Fuel particle dimension and arrangement were 
measured for individual clumps; fuel distribution, quantity and moisture content were 
measured using line transect methods. Bare ground between clumps was also measured. 
Wind speed and direction, air temperature and relative humidity were measured at 10- 
min intervals. Wind speed ranged over 1.11 - 10 m s _1 and FMC over 12-31%. Fire spread 
was measured using metal markers placed near the flame front at intervals of 1-4 mins 
and later surveyed. Fires were allowed to spread until they self-extinguished. The range 
of ROS was 0-1.53 m s^ 1 . Data gathered were analysed using multiple linear regression 
techniques. 



Burrows et al.l (Il99ll ) found that above a threshold wind speed zone (3. 33-4.72 m s *), in 
which flames are tilted sufficiently to bridge the gap between hummocks ( Bradstock and Gilll . 



19931 : iGill et all Il995l ). the ROS varied with the square of the wind speed (R 2 = 0.85). 



Below the threshold wind speed zone, which depends on the percent cover of fuel (ratio of 
percentage of area covered by hummocks to bare ground), the fire does not spread. The 
higher the percent cover, the lower the threshold wind speed required. A lesser, negative, 
linear correlation was determined with FMC. Percent cover and air temperature were also 
found to influence the ROS but much less than either wind or FMC. Fuel load and other 
fuel characteristics were found not to be important. 



Canadian Forest Fire Behaviour Prediction (CFBP) System (1992) 

The Canadian Forest Fire Behaviour Prediction (CFBP) System is a comp onent of the 



Canadian Forest Fire Danger Rating System (CFFDRS) (IStocks et al.l . Il99ll ). which also 
incorporates the Canadian Forest Fire Weather Index (CFWI) System. The CFFDRS is 
the result of continuing research into forest fire behaviour since the mid- 1920s and has 
undergone several incarnations in that time. The current CFFDRS system came into being 
in the late 1960s in the form of a modular structure. The first major component to be 
completed was the CFWI in 1971, which provided a relative measure of fuel moisture and 
fire behaviour potential for a standard fuel type, and has been revised several times since 
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its introduction (IVan Wagnerl . 119871 ). Whil e there have been sev eral interim editions of the 
CFBP, the first of which appeared in 1984 (ILawson et allll985l). it was not until 1992 that 
a final version of the prediction syst em was released ([Forestry Canada Fire Danger Group . 
19921 ; iTaylor and Alexander! . 120061 ) and thus is covered in this review. 



The CFBP system, following on from the long-established Canadian approach to studying 
wildland fire, is based on the combined observations of nearly 500 experimental, prescribed 
and wild fires in 16 discrete fuel types covering 5 major groups: coniferous, deciduous, 
mixed wood, slash and grass fuels. The experimental work on which the system is based 
was conducted by individual researchers working in specific fuel types and locales across 
the country using a variety of methods and published in a variety of places (initially 
in cluding the 1966 work of McArthur in Australian grasslands, later r epl aced by the data 
of Cheney et al. ( 1993 ) ( Forestry Canada Fire Danger Group . 1992h ) . Alexander et al. 
(Il99ll ) provides an overview of the methods used since the 1960s to obtain the dataset 
from the CFBP was derived, but which, due to a number of factors (including technological 
improvements) evolved over the years. The result is a system constructed by a small group 
of dedicated researchers over a period of 20 years that has broad applicability to a wide 
range of fuels and climates. 



Experimental burn plots varied in size from 0.1 ha up to 3 ha (I Alexander et all Il99ll ). 
with the majority being less than 1 ha. Ignition methods included both point ignitions 
as well as line-ignitions. Wind speed unaffected by the fire was measured at 10-m in the 
open (or converted a 10-m in the open equivalent). Experiments were usually conducted 
in the late afternoon in order to attain maximum burning conditions for the day. ROS was 
normally measured by visual observations of fire passage over predetermined distances. 
For point ignition experiments, metal tags were placed at the head and flanks of the fire 
and surveyed afterwards. 

The final version of the CFBP system works in conjunction with the CFWI system to 
determine an Initial Spread Index (ISI) for the standard fuel type (pine forest) and based 
solely on fine FMC and wind speed. The functions chosen for the effect of wind speed 
and fine FMC on the ISI are exponential (exponent 0.05039) for wind, and a complicated 
mix of exponential and power law (exponents -0.1386 and 5.31 respectively) for FMC 
(IVan Wagnerl . 119871 ). No quantification of performance of these functions is given. 



To predict ROS, the ISI is modified by a Build-up Index (BUI), which is a fuel-specific 
fuel consumption factor that includes fuel moisture. Predicted ROS is the headfire ROS 
on level terrain under e quilibrium conditions, thereby implici t ly inc luding effects of ac- 
celeration and crowni ng ( Forestry Canada Fire Danger GrourJ. 19921). The eff ect of slope 
(IVan Wagnerl . Il977bl ) and crown fire transition effects (IVan Wagnerl . Il977al ) then mod- 
i fy the basic ROS. Recent work of the International Crown Fire Modelling Experiment 
(jStocks et all 120041 ) has investigated the behaviour of fully-developed crown fires (which 
is not covered in this review as it is outside the scope of surface fire spread). 



Button (1995) 



Marsden-Smedley and Catchpold (Il995bl ) presented a model for the prediction of ROS 
and flame height of fires in Tasmanian buttongrass m oorlands, described as largely tree- 
less communities dominated by sedges and low heaths (IMarsden-Smedley and Catchpold . 
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1995al ). The behaviour of 64 fires (of which 44 were experimental fires, 4 test fires, 11 
fuel reduction fires and 5 wildfires) at 12 sites was measured. Experimental burns were 
conducted on blocks of either 0.25 or 1.0 ha with ignition line lengths of 50 or 100 m 
respectively under a limited range of weather conditions. ROS was measured by either 
using metal tags thrown at different times or by timing the passage of flames past pre- 
measured locations. For experimental fires, wind speed and direction, temperature and 
relative humidity were measured at 10 m, and wind speed only at 1.7 m above ground 
level, all averaged over 1-3 min periods. Meteorological data for non-experimental fires 
were collected using handheld sensors at 1.7 m. Data ranged from 0.19 - 10 m s -1 for 
wind speed and 8.2-96% for FMC, and - 0.92 m s" 1 for ROS. 



Marsden-Smedley and Catchpold (jl995bl ) found surface wind speed, dead FMC and fuel 
age (time since last fire) to be the key variables affecting ROS, with wind being the domi- 
nant factor. Age and FMC each accounted for 15 to 20% of the observed variation in ROS. 
A power law with an exponent of 1.312 was used to describe the effect of wind, whereas 
both the FMC and fuel age were modelled as exponential functions (FMC decreasing, age 
increasing to a maximum at about 40 years). Rates of spread of the back and flank of 
the fires were found to be approximately 10% and 40% of the head ROS, respectively. 



CALM Mallee (1997) 



McCawl (119971 ) conducted a large-scale field experiment in Eucalyptus tetragona mallee- 
heath community in south-west Western Australia. Shrubs < 1.0 m tall comprised more 
than half the plant species present. Burn plots 200 m x 200 m were established in 20- 
year-old fuel in flat terrain. A semi-permanent meteorological site was set up 500 m from 
the experimental plots recording 30 min averages of temperature and relative humidity at 
1.5 m and wind speed and direction at 2 m. During each experiment, mean wind speed 
and direction at a location up to 250 m upwind of the plot were measured at heights of 2 
m and 10 m at 30 s intervals. FMC was measured using 5 samples of four fuel components 
(3 dead and 1 live) collected post-fire within 30 min of ignition. Wind speed at 10 m in 
open ranged from 1.5-6.9 m s _1 , FMC 4-32%. Experimental fires were ignited using a 
vehicle-mounted flame thrower to establish a line perpendicular to the prevailing wind 
up to 200 m long. Fire spread was measured using buried electronic timers (placed on 
a 24-point grid) equipped with a fusible link that melted on exposure to flames. ROS 
ranged 0.13-0.68 m s -1 . 

Isopleths representing the position of the fire front at successive time intervals were fitted 
to the grid of timer data for each plot using a contouring routine based on a distance- 
weighted least squares algorithm. ROS up to 0.67 m/s and fireline intensities up to 14 
MW/m were recorded. Fires were found to spread freely when the FMC of the dead 
shallow litter layer beneath the low shrubs was < 8%. Forward ROS was modelled as a 
function of the wind speed in the open at 2 m and FMC of the deep litter layer. These 
accounted for 84% of the variation in ROS. A power function (exponent 1.05) and an 
exponential (coefficient -0.11) w e re cho sen to describe wind (measured at 2 m) and FMC 



influences respectively (IMcCawl . 119971 . page 142). Good agreement between the model 
and observations of rate of spread of a limited number of prescribed and wild fires (up to 
ROS = 1.1 m/s), although observed ROS of a wildfire burning under extreme fire danger 
conditions was over-predicted by 30%. 
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CSIRO Grass (1997) 



The CSIRO Grassland Fire Spread Meter (ICSIROl. 119970 is a c ardbo ard circular slide 
rule that encapsulates the algorithms developed by I Cheney et al.l (119981 ) for fire spread in 
natural, grazed and eaten-out grassland pastures. These algorithms are based primarily 
on the results of experiments conducted in annual grasses of the Northern Territory with 
the aim of determining the relative importance of fuel c haracteristics on rat e of forward 
spread of large unconstrained fires, particularly fuel load ( Chenev et al. . 19931), augme nted 
by large experimental fires conducted in open woodland (jCheney and Gouldl . Il995l ) and 
detailed observations of 20 wildfires. 121 experimental fires were carried out on a flood 



plain in a range of fuel treatme nts under a variety of weather conditions (j Cheney et al. 



19931 ; ICheney and Gouldl . 119951 ) in prepared blocks ranging in size from 100 x 100 m to 
200 x 300 m. These fires were predominantly lit from lines ranging in length from 30 
to 175 m, although there were also a number of point ignitions, and allowed to burn 
freely. The range of fuel treatments included mowing and removing cuttings, mowing and 
retaining cuttings, or leaving the grass in its natural state. Two distinct grass species 
(Eriachne burkittii and Themeda australis), of different height, bulk density and fineness 
were present. 

Fuel characteristics (height, load, bulk density, etc.) were measured on four transects 
through each plot approximately every 25 m. In addition to remote standard 10 m and 
2 m meteorological stations, the wind speed at 2 m was measured at the corner of each 
plot and averaged for each ROS interval. FMC samples were taken before and after each 
fire. ROS and flame depth were measured from a series of rectified time-stamped oblique 
aerial photographs of each fire. Wind speed ranged from 2.9 - 7.1 m s _1 , FMC 2.7-12.1%, 
and ROS 0.29-2.07 m s -1 . 



Cheney and Gouldl (119951 ) found the growth of the fires to be related to wind speed and 
the width of the head fire normal to the wind direction. They found that the width of 
the fire required to achieve the potential quasi-steady ROS for the prevailing conditions 
increased with increasing wind speed, and the time to reach this quasi-steady ROS was 
highly variable. ROS was found to depend on the initial growth of the fire, the pasture type 
(natural, grazed or eaten-out), wind speed and live and dead FMC. Utilising the notion of 
potenti al quasi-steady ROS and a minimum threshold wind speed for continuous forward 
spread, ICheney et al.l (119981 ) developed a model of fire spread assuming a width necessary 
to reach the potential ROS. This model uses wind speed, dead FMC and degree of curing 
to predict the potential (i.e. unrestricted) ROS for the prevailing conditions. Above a 
threshold of 5 km h _1 the ROS is assumed to have a power function (with an exponent 
less than 1 (0.844 )) relation with the w i nd sp eed. This wind speed function is similar to 
that proposed by iThomas and Pickardl (Il96ll ). in which a power function with exponent 
of just less than 1 was found. Less than the threshold, the ROS is linear with wind speed 
and dominated by dead FMC. 



Heath (1998) 



A coop erative research effort from a number of Australasian organisations, Heath (ICatchpole et al. 



1998al ) utilises observations of 133 fires (comprising a mix of experimental (95), prescribed 
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(22) and wild (16) fires) conducted i n mixed heathland (heath and sh r ub) fue ls. This in 



eludes 48 experiments conducted by iMarsden-Smedley and Catchpold (jl995bl ) in button- 



grass. Only experimental and prescribed fires were used in model development; wildfire 
observations were used for validation. 

In mixed heathland (comprising heath, scrub and gorse in New Zealand and mixed species 
including Banksia, Hakea and Allocasuarina in Australia), fuel age ranged from 5-25 years. 
Fires were lit as lines of unstated length on slopes < 5°. Due to the disparate nature 
of the researchers involved, methods for measuring variables varied from experiment to 
experiment. Wind speed was generally measured by handheld anemometry at 2 m at 
20 s intervals and averaged over the life of the fire. Wind speed ranged from 0.11-10.1 
m s -1 and ROS 0.01-1.00 m s _1 . Fuel load does not appear to have been measured 
but fuel height was. FMC was measured in some cases and modelled in others using 
pre-established functions based on air temperature and relative humidity. 

Wind speed was found to account for 53% of the variation in ROS. Aerial dead fuel (i.e. 
those fuels not in contact with the ground) FMC was found not to be significant. Fuel 
height was highly significant and with wind accounted for 70% of the variation in ROS. 
A power function of wind speed (exponent 1.21) was used to describe this variation. A 
power function was also used for fuel height (exponent 0.54). 

The model was found to perform reasonably well for the selection of wildfires, considering 
the paucity of available data and necessary assumptions about the involved fuel charac- 
teristics (fuel height, moisture etc.) but could be improved with more variables. The 
wind power function does fail for zero wind but was found to better fit the data than an 
exponential growth function. 



PortShrub (2001) 



Fernandesl (120011 . Il998l ) presented a model developed from field experiments and obser- 



vations of prescribed bur ns conduct e d in f our different types of shrub in flat terrain of 



Portugal. He found that iRothermell (119721 ) did not predict observed ROS well. 29 fires 
were conducted on flat (<3°slope) in gorse, low heath, tall heath and tall heath/tree mix. 
Fine aerial live and dead FMC was sampled prior to each burn. Meteorological variables 
(wind speed, air temperature, relative humidity) were measured at 2 m in the open us- 
ing either a fixed weather station placed near the burn plot or upwind with handheld 
instruments. Fires were lit as lines of length 10 m in experimental fires and 100 m in 
prescribed burns. ROS was measured by recording time of arrival of the head fire at 
reference locations. Wind speed ranged 0.28-7.5 m s" 1 , FMC 10-40% and ROS 0.01-0.33 
rn s _1 . 

ROS was significantly correlated with wind speed (1% level) and less so with RH, temper- 
ature, and aerial dead FMC (5%). Other fuel characteristics were also found to affect ROS 
but were strongly intercorrelated and thus could not be separated, however, preference 
was given to fuel height. The initial model found a power law (exponent 1.034) for wind 
speed. However, as the model predicted no ROS in zero wind, an exponential function 
(coefficient 0.092) was subsequently incorporated. The final model, with an exponential 
decay function for dead FMC (coefficient -0.067) and power function (exponent 0.932) for 
fuel height, improved the overall performance of the model (R 2 = 0.91). The model was 
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also found to predict we ll the data sets of other authors and be i n close agreement wit h 
other field studi e s (e.g. iMarsden-Smedley and Catchpold (jl995bl ); I Cheney et al.l (119981 ); 
Catchpole eiTa!! ril998a[ )). 



CALM Jarrah I (1999) 



Burrowsl fll999al . Il994h conducted a series of 144 laboratory experiments (54 wind-driven, 



6 no wind, 13 backing, 34 with slope, 15 point ignition) using fallen leaves and twigs (<6 
mm) placed on a 4 m long by 2 m wide table set in a large shed. Wind was supplied 
by four domestic fans calibrated to give a desired wind speed over the fuel bed. FMC 
was varied by uncontrolled ambient conditions and wetting prior to burning. It ranged 
from 3 to 14%. Fires were lit along the 2 m upwind edge using cotton wick soaked in 
methylated spirits and allowed to burn for 50 cm before measurements commenced. ROS 
was measured by recording time taken to reach end of fuel bed. Wind was varied from 
0.0 to 2.1 m s" 1 with mean 1.06 m s" 1 . ROS ranged 0.002-0.075 m s" 1 . 

In wind-driven fires, no relationship between fuel load and forward ROS was found. Most 
variation in ROS was due to wind speed (correlation coefficient 0.94). ROS was negatively 
related to FMC (correlation coefficient -0.31). Backing ROS was found to be directly 
related to fuel load. 

At wind speeds < 0.83 m s _1 , ROS was relatively insensitive to wind. Above this value, 
ROS was found to vary linearly with wind speed. However, a power function (exponent 
2.22) was used to model wind speed effect on ROS. An inverse linear function was used 
for FMC. This model was found to underspecify ROS > 3.33 m s _1 with an error variance 
that increased with ROS. 



CALM Jarrah II (1999) 



Burrowsl (jl999bl . Il994l ) studied four series of fire behaviour data obtained from field ex- 
periments and fuel reduction burns on flat to gently sloping terrain in Jarrah (Eucalyptus 
marginata) forest in south-west Western Australia to test Jarrah I and other models for 
forest fire spread. Fuel was characterised by a layer of dead leaves, twigs, bark and floral 
parts on the forest floor with low (<0.5 m, 30% cover) understorey of live and suspended 
vegetation. Plots were 100 m wide x 200 m long. 56 of 66 total plots were lit from lines 
of 50-100 m length, the remainder being point ignitions. 

Historical (pre-fire) weather data (including rainfall, temperature, relative humidity, wind 
speed and direction at 2 hourly intervals) were obtained from nearby permanent weather 
stations. During each experiment a portable weather station approximately 50 m from 
the fire recorded wind speed at 1.5 m and 10 m and temperature and relative humidity 
at 1.5 m at 5 minute averages. FMC was measured at the time of ignition. Wind speed 
at 10 m in the open ranged 0.72-3.33 m s _1 and FMC 3-18.6%. 

ROS was measured by recording the time of arrival at a grid of predetermined locations, 
along with other fire characteristics, after first allowing the fire to spread 20 - 40 m (~ 
15 min) in order for it to attain a quasi-steady ROS for the prevailing conditions. The 
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position of the flames in relation to the grid was mapped at 5 min intervals. ROS ranged 
0.003-0.28 m s _1 . 



Unlike the laboratory findings ( Burrows . 1999al ). Burrows here found a non- linear relation 



between wind speed and ROS. A power function (exponent 2.674) was selected. FMC was 
determined to also be a power function (exponent -1.495). Like the laboratory findings, 
fuel load was not found to correlate with ROS. The model was found to underpredict 
ROS of large wildfires burning under severe conditions. 



Gorse (2002) 



Baeza et al.l (120021 ) conducted field experiments during spring and autumn in gorse shrub- 
lands of eastern Spain with the aim of developing a prescribed burning guide. Fuels were 
3, 9 and 12 years old and were replicated 3 times resulting in a total of 9 fires. Plots were 
33 m x 33 m and were burnt under low (< 1.39 m s _1 at 2 m) wind, utilising headfire 
spread for the 3-year-old fuel and backing fires for the other two age classes. Meteoro- 
logical data was recorded at 2 m at 15 min intervals. Fuel characteristics were recorded 
along 5 parallel transects 5 m in length. FMC was measured from 10 samples of the most 
abundant species collected prior to ignition and ranged from 22-85%, presumably includ- 
ing both live and dead fuelsj. Ignition technique is not specified. ROS was measured by 
recording the time to travel a fixed distance within the plot and ranged 0.004-0.039 m 



It was found that FMC was the dominant factor affecting ROS in a linear manner (coeffi- 
cient 0.487). The combination of heading and backing propagation negated any consistent 
effect of wind speed on ROS. 



PortPinas (2002) 



Fernandes et al.l (120021 ) developed a model for the behaviour of fires in maritime pine 
(Pinus pinaster stands in northern Portugal under a range of fire weather conditions 
that occur outside the wildfire season for the purpose of improving the understanding of 
prescribed fire for hazard reduction. Six study sites in mountainous terrain with forests 
founded by plantation or regeneration following fire events and aged 14 to 41 years were 
established. Fuel complexes were dominated by litter, shrubs or non-woody understorey 
(e.g. grass) types. Extensive destructive and non-destructive sampling to quantify the 
fuels was undertaken along transects in each experimental plot. Four strata of fine fuel 
layers were defined: shrubs, herbs and ferns, surface litter and upper duff. Experimental 
plots were square, 10-15 m wide, and defined by 0.3 to 1.2 wide control strips assisted by 
a hose line during burning. 

Wind speed was measured continuous at 1.7 m above ground approximately 10 m from 
each experimental plot. Three composite fuel moisture samples (one litter, one duff and 
one live) were sampled at random locations prior to ignition. 



5 It should be noted that the authors dri ed their fuels at 80° for 24 hours which is much less than the 
generally accepted 104°C for 24 hours (e.g. IChenev et all (|l993h ). perhaps resulting in lower than actual 



FMC values-see discussion on Measurement issues. 
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94 experimental fires for fire behaviour studies were conducted when slope and wind 
direction were aligned within 20°. Line ignition occurred 2 m from the windward edge to 
allow both forward and backing spread observations. Fire behaviour measurements used 
1.5-m-high poles located at regular distances along the plot reference points. ROS 

was determined by recording the time at which the base of the fire front reached each 
pole. Flame height and flame angle were estimated visually and used to calculate flame 
length. Wind speed ranged from 0.3 - 6.4 m s _1 , surface dead FMC ranged 8 - 56%, air 
temperature 2 - 22°C and relative humidity 26 -96%. ROS ranged 0.004 - 0.231 m s _1 . 



Fernandes et al.l (120021 ) found that three existing models underestimated ROS with sig- 
nificant differences between predicted and observed values, as much as 8-fold in one case. 
Undertaken non-linear least=squares analysis, they found that slope and wind speed were 
the most significant variables with dead FMC in a less significant role. A power law func- 
tion with wind speed only (exponent 0.803) explained 45% of the variation in ROS. If 
wind speeds less than 0.83 m s _1 were excluded, the correlation coefficient increased to 
0.996. Slope alone explained 30% of the variation. The final model selected for litter- 
shrub fuels (the general case for maritime pine stands) involved wind speed (power law, 
exponent 0.868) dead surface FMC (exponential, coefficient -0.035), slope and understorey 
fuel height. Fuel height was selected as could be considered a surrogate for the overall 
fuel complex structure effect on ROS. The model was then adapted through changes in 
constants to predict ROS in litter and non-woody understorey complexes. No assessment 
of the performance of this model was reported. 



Maquis (2003) 



Bilgili and Saglaml (120031 ) conducted a series of 25 field experiments in open, level shrub- 
land of maquis fuel in southwestern Turkey. Average height of the fuel was 0.53 m and 
fires were conducted under a range of wind and fuel conditions. Each fire plot was 20 m 
wide by 30 m long. A meteorological station recorded air temperature, relative humidity, 
wind speed and precipitation at 1.8 m daily. Fuel characteristics were measured from 
random destructive sampling prior to the experiment series. Live and dead FMC was 
sampled immediately prior to ignition. During each fire, wind speed, temperature and 
relative humidity at 1.8 m were recorded at 1 min intervals using the automatic meteo- 
rological station. These were averaged over the period of fire spread. Wind speed ranged 
0.02-0.25 m s~\ FMC 15.3-27.7%. 

Fires were lit with a drip torch along the upwind (20 m long) edge to quickly establish a line 
fire and were allowed to propagate with the wind across the plot. ROS was measured by 
recording time of arrival at a series of predetermined locations and ranged 0.01-0.15 m s _1 . 
ROS was strongly correlated with wind speed; a linear function explained 71% of observed 
variation. FMC was found not to have any significant effect on ROS, attributed to the 
narrow range studied. The final model used a linear function of wind speed (coefficient 
0.495) and total fuel load, with an R 2 = 0.845. 
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Quasi-empirical models 



Where the data gathered from experimental observation is analysed using a physical 
framework for the functional relationships between dependent and independent variables, 
a quasi-empirical model results. The degree to which the physical framework controls the 
structure of the model can vary but the nature of the model is essentially based upon the 
observed data (which differentiates it from quasi-physical models which use data solely 
for parameterisation). Table [5] summarises the quasi-empirical models discussed below. 



TRW (1991) 



Wolff et al.l (119911 ) presented the results of laboratory experiments conducted in a purpose- 
built wind tunnel 1.1 m wide by 7 m long with a moveable ceiling. The fuel layer was 
vertical match splints (1.3-4.4 mm in diameter) set in a ceramic substrate. Wind speed 
varied from 0-4.7 m s _1 , ROS rang ed from 0-0.0 7 m s -1 . The results confirmed the 



theoretical treatment conducted by ICarrier et al.l (1l99ll ). in which it was hypothesised 



that the dominant heat transfer mechanism in such a set-up would be a mix of convection 
and diffusion (i.e. 'confusion') heating that would result in a relationship in which the ROS 
would vary as the square root of the wind speed normalised by the fuel load. If radiation 
was the predominant preheating mechanism, it was hypothesised that the variation would 
be as the power of 1.5 rather than 0.5. 

Wolff et al. found that not only did the width of the fuel bed play an important part in 
determining the ROS but also the total width of the wind tunnel itself. The narrower 
the fuel bed, and the facility, the slower the ROS. It was suggested that a narrower fuel 
bed forced air away from the fuel bed due to drag considerations in the fuel. A series of 
experiments with tapering fuel beds and working section confirmed this. If the fuel bed 
and working section was too narrow, ROS ceased. 



NBRU (1993) 



Beerl ( 1l993bl . Il99ll ) investigated the interaction of wind and fire spread utilising a series 
of 18 small-scale wind tunnel (length 40 cm, height 16 cm) experiments using a single 
row of match splints in wind ranging from 0.0 to 9 m s -1 . ROS ranged from 0.004-0.38 m 
s _1 . Rather than a single continuous function to describe the relationship between wind 
speed and ROS, Beer put forward the hypothesis that there exists a critical characteristic 
(threshold) wind speed that affects ROS with different wind speed functions above and 
below this value. Below the threshold Beer found a normalised (by the threshold wind 
speed) power function (exponent 0.5). Above the threshold, Beer found a normalised 
power function (exponent 3.0). Beer postulates that the choice of the value is related to 
the wind speed at which the wind shear is strong enough to generate flame billows and 
that this value corresponds to a mid-flame wind speed of 2.5 m s" 1 . Above this value it is 
thought that the flames remain within the fuel bed rather than above it. Beer attempted 
to fit this model to observations of grassfire behaviour but could not. 



Beerl (Il993al ) further explored the effects of wind on fire spread through simplified (match 



splints) fuel. His extension of a simple geometric model of fire spread in no wind to include 



15 



wind (based on geometry of wind-tilted flame and distance between fuel elements), in 
which the ROS-wind function is a complicated solution to a set of equations to determine 
the critical time for flame immersion of adjacent fuel elements, did not perform well. 
This was attributed to assumptions about the characteristics of the flame and a constant 
ignition temperature. Beer concludes that a single simple power law or exponential is 
unlikely to be a correct mathematical description for the ROS-wind speed relation. 



USFS (1998) 



Catchpole et al.l (jl998bl ) conducted an extensive series of environmentally-controlled wind 
tunnel experiments and used the results, in conjunction with energy transfer considera- 
tions, to develop a spread model, USFS (United States Forest Service). 357 experimental 
fires were carried out on a fuel bed 8 m long by 1 m wide in a 12-m long wind tunnel 
of 3 m square cross section. Four fuels with different surface-area-to-volume ratios (two 
sizes of poplar excelsior, ponderosa pine needles and ponderosa pine sticks) were chosen to 
be reasonable approximations to natural fuel layers. Temperature and relative humidity 
were controlled to produce a range of FMC, 2% to 33% (although the majority of fires 
were carried out at ambient values of 27°C and 20% RH giving an FMC range of 5-9%). 
Wind speed above the tunnel's boundary layer ranged from 0.0 to 3.1 m s _1 . Rate of 
spread was measured at 0.5 m intervals using photovoltaic diodes placed 25 mm above 
the fuel bed to record the time of arrival of the flame front. 



Utilising the conservation of energy model of lFransdenl (119711 ). Fransden's (119731 ) effective 
heating number, a propagating flux model that is linear in packing ratio, and an expo- 
nential decay function for FMC, the author s built a model of fire spread very similar in 
its construction to that of iRothermell (119721 ) except that they used the heat of ignition of 
a unit mass of fuel (which comprises the heat of pyrolysis and heat of dessication) rather 
than the heat of pre-ignition as used by Rothermel. A power function for wind was then 
fitted to the data and an exponent of 0.91 determined. Although a cubic polynomial 
function was found to better fit the data, the authors chose the power function as it was 
more consistent with data from wildfire observations. 



Coimbra (2002) 



Viegasl (120021 ) presents a quasi-empirical model of fire spread that utilises the geometry of 
the fire perimeter to de termine the fo rward spread rate. The m ain conceit of this notion, 
previously proposed in IViegad (119981 ) and IViegas et al.l (119981 ). is that a line fire lit at 
an angle to a slope or wind gradient undergoes a translation and rotation of the fireline 
in order to spread with the maximum rate in the direction of the gradient. Extensive 
laboratory experimentation utilising a double-axis tiltable fuel bed (1.6 m x 1.6 m) for 
a range of forward/back and left/right slopes was used to develop the model. The fuel 
was Pinus pinaster needles with an FMC determined by ambient conditions (ranging 
10% - 15%). 23 experimental fires we r e con ducted, with 10 fires of varying slope and 
13 fires of varying inclination. IViegasl (120021 ) found a maximal rotation velocity at an 
inclination angle of 60° but was unable to convert this to a forward ROS. However, 
Viegas does develop a fire perimeter propagation algorithm in which the perimeter is 
treated as a continuous entity that will endeavor to align itself with the gradient, through 



16 



this proposed rotation mechanism, to an angle of approximately 60°. The translation 
and rotation hypothesis, however, ignores a basic observation of the evolution of flanking 
spread and instead assumes that spread at non-parallel angles to the slope or wind gradient 
must be driven by a headfire. 



Viegasl (120051 ) attempts to extend these ideas to describe the phenomenon of 'fire blow-up' 
based on the concepts of fire 'feedback effects'. Viegas proposes the existence of a positive 
dynamic feedback between the ROS of a fire and the flow velocity driving the fire such 
that the fire accelerates exponentially. He uses some of the results of experimental fires 
burnt in a "canyon" , a doubly-sloped tray, in no wind and a range of canyon slopes and 
inclinations to parameterise his model and the remainder to test it. Viegas treats all data 
for all slope and inclination combinations as independent and continuous. As a result, 
his model increases ROS exponentially, resulting in extremely rapid acceleration-what he 
describes as blow-up. However, categorised by slope, rather than treated continuously, 
the ROS data actually asymptotes to a reasonable number in each case, which in most 
cases confirms the long-he ld rule of thumb of doubling the flat g round ROS fo r every 10 
degrees increase in slope ( iMcArthurl . 119671 ; IVan Wagnerl . I19881 ) . I Vie gad (120061 ) conducts 
a parameteric study of this model and determines that fires in light and porous fuels are 
more likely to exhibit 'eruptive' behaviour than fires in heavy and compacted fuels. 

Viegas's extrapolation of this model to fatal wildfire incidents is tenuous at best and 
really only proves the widely accepted acceleration up a slope. Other, more robust, 
theo ries of unexpected fir e behaviour resulting in fatalities are probably more applicable 
(e.g. IChenev et all (1200 lh I 



Nelson (2002) 



Nelson! (120021 ) extende d the quasi-empirical wor k of iNelson and Adkinsl (119881 ) utilising 
the laboratory data of I Weise and Bigingi (119971 ) to build a trigonometric model of fire 
spread that combines wind and slope effects into a single combined 'effective' wind speed. 



T he INelson and Adkinsl (119881 ) model utilised the dimensional analysis of fire behaviour 
of Byram ( 19661 ). where three dimensionally homogeneous (i.e. dimensionless) relations 
were derived: 1) the square root of the Froude number, 2) a buoyancy number relating 
convective heat output to rate of buoyancy production, and 3) t he ra tio of combustion 
time to time characteristic of flame dynamics. INelson and Adkinsl ( 119881 ) then used spread 
observations from 59 experimental fires (a total of 44 lab and 21 field, some deleted) and 
mixed and matched the dimensionless relations until they found a combination that gave 
a reasonable correlation. They derived a dimensionless form of ROS and wind speed, 
which, when fitted to the data and converted back to dimensions, gave a power law 
relation between wind speed and ROS (exponential 1.51). As the maximum wind speed 
used to obtain the data was 3.66 m s -1 and maximum observed ROS was 0.271 m s -1 , 



Nelson and Adkinsl (119881 ) acknowledged the need for higher wind speed experiments. 
ROS was also found to be a function of fuel load (power function, exponent 0.25) and 
residence time (inversely proportional). FMC was considered to be accounted for in the 
estimate of fuel load and residence time. 



Rather than the traditional approach used bv lMcAlpine et al.l (119911) where the equivalent 
wind speed for a slope-only ROS was determined, INelson! (120021 ) used the concept of 
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vertical buoyant velocity and the slope angle component of this wind to construct a 
slope-induced wind which was then added yectorally to the ambient wind across slope. 
Nelson extended the dimensional analysis of iNelson and Adkinsl (119881 ) to then determine 
ROS. The resultant equations, which do not apply to flanking or backing fires due to the 
assumption about connective heating throu gh the Froude number, were then compared 
against the data of IWeise and Biging) (119971 ). gathered from 65 experiments in a portable 
tilting wind tunnel using vertical paper birch sticks as the fuel bed in a variety of wind 
and slope configurations, ranging from 0.0 to 1.1 m s _1 and -30°to +30°. The effective 
wind speed was found to correlate linearly with ROS. 



Discussion 



Wind speed function 

As stated earlier, one method of comparing the structure of each of the above models is 
to examine the form of the functional relationship between ROS and wind speed chosen 
by the authors. Table [3] summarises the models discussed and the form of the wind speed 
function chosen. Also listed are the experimental bounds of the wind speed and ROS. 

Only three of the models for which the wind function is given are not power functions, 
PortShrub and CFBP are exponential and Maquis is linear. Of the remaining models, 
3 models have exponents less than one: TRW, CSIRO Grass and USFS. The remaining 
wind functions all result in non-linear increases (Figure d]) in the R OS that will result in 
a speed greater than the wind speed driving it, which is unphysical (Beer, 199 lh . (While 



this is also the case also for CFBP as illustrated here for wind speeds > 15 m s _1 , the ROS 
function is further modified in the CFFDRS system by fuel-specific functions which can 
reduce the predicted ROS below the wind speed.) The reason for this choice of function 
appears to be the desire by the modellers to fit data at low wind speeds (including zero). 
Many of the models had ranges of wind speed that were fairly low (< 3 m s" 1 ). 

The few models that were based on large ranges of wind speed in field experiments (with 
the exception of CALM Jarrah II) tended to result in power functions with exponents less 
than one. CSIRO Grass is the highest power function less t han one and is very simila r in 



form to the linear function of Maquis over the given range. iFendell and Wolfll (120011 ). in 



their brief review of the topic including a n umber of older models, found that the w ind 



power function exponents ranged from 0.42 (IThomasl . 119671 ) to 2.67 (jBurrowsl . Il999bl ). 



There seem to be two key factors in the choice of functional relationship used to describe 
fire spread and wind speed. The first is the need to fit the function through the origin. In 
many cases, particularly laboratory experiments, zero wind speed is taken as the default 
state and thus any continuous function must not only fit the data of non-zero wind, but 
also of zero wind. This is discussed in greater detail below. The second is that for the 
most part, the range of wind speeds studied (again particularly in the laboratory) is very 
small. As can be seen in Figure [TJ any function, be it cubic or very shallowly linear, 
performs rather similarly at low wind speeds (<1.5 m s _1 ). 

It is interesting to note that in the full range of functions presented here, the nearly 
median wind function in Figure [T] (i.e. Heath) is the result of the combination of multiple 
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datasets, experimental methods and authors, perhap s resulti n g in a midd l e gro und of 
approaches. Many physical models of fire spread (e.g. iGrishinl (119841 ); iLinnl (119971 )) have 
observed linear functional relationships between wind speed and ROS, suggesting that 
power law functions with exponents close to unity may have a more fundamental basis. 



In their validation of the performance in Mediterrane an shrub fuels of seven wildland fire 
spread models, including the CFBP and Rothermel, ISauvagnargues-Lesage et al.l (120011 ) 
found that a model's performance is not related to the model's complexity and that even 
the most simple model (in this particular case a local fire officer's rule of thumb based on 
a linear discount of the wind speed) performed as well as more complicated models such 
as CFBP. 



Threshold wind speed 



One important aspect differentiating the various choices of wind function, is the ideal 
of a continuous funct ion that includes zero wind speed. It has been noted previously 
(IBurrows et al.l . Il99ll ) that fires in discontinuous fuels such as spinifex have a minimum 
wi nd speed requi r ed bef ore forward spread is achieved. This notion was extended further 
by I Cheney et al.l (119981 ) to define a threshold wind speed at which fires spread forward 
consistently. The argument is that fires burning in low winds in the open respond to 
eddy fluctuations in the wind flow (resulting in near circular perimeter spread after a long 
period) and do not spread in a continuous consistent manner until the wind speed exceeds 
a certain threshold. Above this threshold, the fire spreads forward in a manner directly 
related to the wind speed. 



The choice of threshold value is dependent then upon the method of measuring the wind 
speed (location, height, period, etc) and the fue l type in whic h the fi re is burning (taller 
fuels reduce the wind speed reaching the fire). IChenev et al.l (|l998T) ch ose a 1.39 m s -1 
open wind speed threshold for open grasslands. iFernandes et al.l (120021 ) found that wind 



speed explained more variation in ROS when wind speeds below 0.83 m s" 1 (at 1.7 m) 
were excluded, suggesting that at low wind speeds factors other than wind play a more 
significant role in determining ROS. 



Fuel moisture content function 



Another method of comparing the various empirical and quasi-empirical models is that 
of fuel moisture content function. Not all the models discussed here addressed the rela- 
tion between fuel moisture content and rate of spread, and so the discussion here is not 
comprehensive. Figure [2] shows the various functions for those models that include fuel 
moisture content. 

As with the wind function, there is a wide spread of functional forms used to described 
the influence of FMC on ROS, perhaps reflecting modelling approaches, methods or per- 
sonal choice. There appear to be three types of functions representing the fuel moisture 
content/ROS function: weakly linear (e.g. Gorse (normalised) or CALM Spinifex (nor- 
malised)) in which the FMC plays a minor role in determining the ROS, strongly expo- 
nential (e.g. CALM Jarrah II and USFS) in which FMC plays a strong role until very low 
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FMC values, and strongly linear or weakly exponential (which might be approximated to 
linear) in which the role of FMC is spread over a large range of values. The majority of 
models discussed here fall into the latter group. 

The weakest of the linear models (Gorse (normalised) is characterised by few experiments 
with a limited range of FMC values, which raises the issue of how many sample points are 
needed to properly inform functional choice and model validation-one could argue that 9 
fires is simply not enough given the range of the uncontrolled variables. 

The weakly exponential group of models, which includes strictly linear models (e.g. 
CALM-Spinifex), appears to be the most robust in terms of range of FMC values and ex- 
perimentation. It is interesting to note the similarity between the CALM Mallee and the 
CSIRO Grass models. The large difference in functionality between the strongly exponen- 
tial and weakly exponential is interesting and may reflect differences in functionality as a 
result of wind function modelling, as all modelling identified wind speed as the primary 
variable and FMC as the secondary. 



Measurement issues 



All empirical science is limited by the ability to measure necessary quantities, to quantify 
the errors in those measurements, and then relate those meas urements to the phenomenon 



under investigation and wildland fire science is no different. ISullivan and Knightl (120011 ) 
discussed the determination of the errors in measuring wind speed under a forest canopy 
some distance from an experimental fire and relating that measurement to measurements 
of fire spread. The issues of where to measure (location, height, in the open, under the 
canopy, etc.), how long to measure (instantaneous, period sampling, average, period of 
average, etc.), and how to correlate measurements with observations are complex and 
necessarily require approximations and simplification in order to be undertaken. 

Similarly, destructive sample of FMC has issues that complicate a seemingly simple quan- 
tity. The time of sampling (morning, afternoon, wetting period, drying period), the gen- 
eral location (in the open or under the canopy, in the sun, in the shade, in between, etc.), 
the specific location (surface litter fuels, profile litter fuels, mid-layer fuels), the species 
of fuel (predominant fuels, non-predominant fuel, live, dead, etc.). Also once the samples 
have been taken there then is the issue of best drying methods for the particular samples 
to ensure a water- free weight, the best method of determining an average value for a plot, 
variance, error, etc. 

Quantifying other factors such as fuel, again seemingly simple quantities rely upon a 
knowledge of the mode of combustion of the fire and which aspects of the fuel most 
influence that combustion and therefore the behaviour of the fire. These include definition 
of fuel strata (which itself depends on the intensity of the fire and which parts of the fuel 
complex will be burning during the fire and thus contribute to the energy released by the 
fire), the structure of the fuel and the size of fuel particles important to fire behaviour 
of the front, flanks and behind the flame zone, the amount of fuel available, the amount 
consumed, the chronology of the consumption of the fuel, the mode of consumption, 
transport of burning fuel (i.e. firebrands), spatial and temporal variation of these fuel and 
fire characteristics, determination of averages and methods of averaging, determination of 
errors, etc. The list could continue. Other factors, such as air temperature and relative 
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humidity, insolation, atmospheric stability, slope, soil type and moisture, have their own 
range of measurement difficulties, and are by no means the only quantities involved in 
quantifying the behaviour of wildland fires. 

Laboratory-based experiments may aim to reduce the variation and control the errors 
in measurement of many of these quantities but are not immune to the difficulties of 
measurement. 



Field versus laboratory experimentation 



Empirical or quasi-empirical modelling of fire behaviour has resulted in significant ad- 
vances in the state of wildfire science and produced effective operational guides for de- 
termining the likely behaviour of wildfires for suppression planning purposes. Unlike 
physical or quasi-physical models of fire behaviour, these systems are simple, utilise read- 
ily available fuel and weather input data, and can be calculated rapidly. However, there 
is a significant difference between those models developed from field experimentation and 
those developed from laboratory experimentation. 

Large-scale field experiments are costly, difficult to organise, and inherently have many of 
the difficulties associated with wildfire observations (e.g. spatial and temporal variation of 
environmental variables, uncontrolled variations, changing frames of reference and bound- 
ary conditions, etc.). Laboratory experiments can be cheap and safe, provide relatively 
repeatable conditions, and can limit th e type and range of variations within variables 
and thus simplify analysis. IVan Wagner! (119711 ) raises the issue of laboratory versus field 
experimentation but avoids any categorical conclusions (perhaps because there are none), 
simply stating some features of wildland fire behaviour are better suited to studying in the 
controlled environment of the laboratory or could not be attempted in the field, and other 
features cannot be suitably replicated anywhere but in large-scale outdoor experiments. 

Correct scaling of laboratory experiments (and field experi ments fo r that matte r) is vital 



to rep licating the conditions expected during a wildfire. iByraml (119661 ) and IWilliams 



■7? 

(jl969l ) conducted dimensional analysis of (stationary) mass fires in order to develop scaling 



laws to conduct scaled model experiments. Both found that scaling across all variables 
presents considerable difficulty and necessitates approximations, particularly in regard to 
atmospheric v ariables , whic h result in impractical lower limits (e. g. model forest fires ~ 
6-16 m across ( Byram . 1966 ). or gravitational acceleration ~ 10 g ( Williams . 1969h ) . As a 
result, it is clear that any scaled experiments must take great care in drawing conclusions 
that are expected to be applicable at scales different from that of the experiments; not 
only may physical and chemical processes behave differently at different scales but the 
phenomena whole may behave differently. 

The key difference between field-based and laboratory-based experimentation, in this 
author's opinion, is the assumptions about the nature of combustion (including heat 
transfer) that are required in order to design a useful small-scale laboratory experiment. 
That is, there is the presumption implicit in any laboratory experiment that there is 
sufficient understanding about the nature of fire such that key variables can be isolated 
and measured without regard to the fire itself. 



One such aspect identified by I Cheney et al.l (119931 ) is the importance of the size and shape 
of the fire in determining resultant fire behaviour. Prior to this work, it was thought that 
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the size of the fire played little part in determining the behaviour of a fire and thus the 
results of small experimental fires could be extrapolated to larger fires burning under 
less mild conditions. Other factors such as the physical structure of the fuel or moisture 
content of live fuels, or other hitherto unconsidered factors, may play less significant but 
important roles in explaining the unaccounted variation in ROS. 

Field experiments on the other hand, by their very nature, are real fires and thus in- 
corporate all the interactions that define wildland fire. This aspect holds considerable 
weight with end users who endow such systems with a confidence that pu rely theoretical 



or laboratory-only-based models do not receive. As iMorvan et al.l (12004 ) concluded, no 
single approach to studying the behaviour of wildland fire will provide a complete solution 
and thus it is important that researchers maintain open and broad paradigm 
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Table 1: Summary of empirical models discussed in this paper 
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Table 2: Summary of quasi-empirical models discussed in this paper 
Model Author Year Country Field/Lab Fuel type No. fires Size (w x 1) 
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Figure 1: Graph of functional relationships between wind speed and rate of forward spread 
used in various empirical and quasi-empirical fire spread models. The relationship is only 
indicative of effect on ROS as the full model in each case may also include effects of other 
variables such as fuel moisture content as well as increases or decreases in wind speed due 
to measurement at different heights. 
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Figure 2: Graph of functional relationships between fuel moisture content and a rate of 
forward spread factor used in various empirical and quasi-empirical fire spread models. 
A number of models have been normalised in order to present them in conjunction with 
other models (i.e. norm.) The ROS factor shows the effect fuel moisture content has on 
the final rate of spread value in each model. 
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